Strong Phase Separation in a Model of Sedimenting Lattices 



> 
O 



o 

OS 



i 

o 



X 



Rangan Lahiri 1 0, Mustansir Barma 1 !*], and Sriram Ramaswamy 2 
1 Department of Theoretical Physics, Tata Institute of Fundamental Research, Homi Bhabha Road, Mumbai 400 005, India; 
2 Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560 012, India 
0\ '■ (February 1, 2008 TO APPEAR IN PHYS REV E, 1 JANUARY 2000) 

G\ 



We study the steady state resulting from instabilities in crystals driven through a dissipative 
medium, for instance, a colloidal crystal which is steadily sedimenting through a viscous fluid. The 
problem involves two coupled fields, the density and the tilt; the latter describes the orientation of 
the mass tensor with respect to the driving field. We map the problem to a 1-d lattice model with 
two coupled species of spins evolving through conserved dynamics. In the steady state of this model 
' each of the two species shows macroscopic phase separation. This phase separation is robust and 

survives at all temperatures or noise levels — hence the term Strong Phase Separation. This sort of 
phase separation can be understood in terms of barriers to remixing which grow with system size 
and result in a logarithmically slow approach to the steady state. In a particular symmetric limit, it 
^ ■ is shown that the condition of detailed balance holds with a Hamiltonian which has infinite-ranged 

interactions, even though the initial model has only local dynamics. The long-ranged character of the 
interactions is responsible for phase separation, and for the fact that it persists at all temperatures. 
Possible experimental tests of the phenomenon are discussed. 
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I , I. INTRODUCTION 

Q ■ A. Background 

,o ■ 

Sedimentation - the settling of heavier particles in a lighter fluid - is a rich source of intriguing physics jjj. The 
steadily sedimenting state arises, of course, from a balance between gravity and viscosity. Viscous damping in this 
^ • noncquilibrium steady state has important consequences: when a given particle is slowed down by the fluid, its 
momentum does not disappear, but produces disturbances in the fluid which affect the motion of other particles 
This makes sedimentation a challenging problem in the statistical physics of driven many-body systems. 

In the general area of noncquilibrium steady states, much recent progress has come by stepping away from the dif- 
ficulties of hydrodynamics and focussing instead on simple driven lattice-gas models In fact, intimate connections 
were discovered by two of the present authors (hereafter LR) ^ between these models and the physics of sedimenting 
Q\ , crystalline suspensions (as well as a closely related problem, a flux-point lattice moving through a superconducting 
slab). The LR model was based on two crucial properties of collective settling discovered by Crowley 0] in his theoret- 
ical and experimental studies of hard spheres sedimenting in a viscous medium: (i) The magnitude of the local settling 
velocity of a region of the crystal depends on its concentration, i.e., on the particle number density in that region, 
and (ii) the direction of the local settling velocity depends on its tilt, that is, the orientation, relative to the applied 
force (gravity) of the principal axes of the local particle distribution. These effects, which also follow from symmetry 
arguments, were incorporated into a natural one-dimensional model for the coupled, stochastic, local spin-exchange 
, dynamics of two sets of Ising variables, {o~i} with states denoted by + and — for the concentration relative to the 
mean, and {r^} with states denoted by / and \ for the tilt, on the sites i of a one-dimensional lattice. Analysis of this 
model leads to several interesting results, some published in H and some new ones which we present here. 



B. Update rules 



Our results will be easier to understand after a quick summary of the update rules of the lattice model, which we 
turn to next. This will also serve to underline the simple nature and potentially wide applicability of the model. It 
is convenient to place the {o~i} and {r^} on two sublattices of our one-dimensional lattice; we label sites on the first 
sublattice by integers, and those on the other by half-integers. A configuration is then a string ri o\ r 3 oi T5 o~sTT_o~i . . ., 
say / + \ — / — / + \ Using the above notation for the states of the two variables, and denoting the rate of 
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an exchange process by W, the probabilities per unit time for the various possible exchanges can be represented 
succinctly by 
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where the first line, for example, represents the rate of H — going to — h in the presence of a downtilt \, and so 
on. The quantities D, E, D' , E' (all positive) and a, b, a', b' are all in principle independent parameters but we will 
argue below that the case of physical interest and relevance to the sedimentation and driven flux-lattice problems 
is sgna = sgna', sgn& = sgnfe', and that the quantity which controls the qualitative behaviour of the model is then 
a = sgn(afe). We find two completely distinct kinds of behaviour, depending on whether a is positive or negative. 
If a < 0, the steady state of the model is a mixture of pluses and minuses, and of uptilts and downtilts, which is 
statistically homogeneous on a coarse-grained level. If a > 0, such a state is unstable with respect to fluctuations 
which drive it to a strongly phase separated state of a type defined and discussed below. We refer to the cases a < 
and a > as the Stable and Unstable LR (SLR and ULR) Models respectively. 



C. Strong Phase Separation: Summary of Results 

The focus of this paper is the study of phase separation phenomena of a new and unusual sort, in the unstable LR 
(ULR) model of sedimenting colloidal crystals described above. Following the appearance of the LR model, the same 
type of phase separation was shown to occur in a three-species permutation-symmetric model on a one-dimensional 
lattice with periodic boundary conditions, with dynamics which may be regarded as a simplified version of that in 
||. A further generalisation which breaks the permutation symmetry between the three species was studied in 
The underlying mechanism of phase separation appears robust and simple enough that it might be worth looking for 
in other systems. Here is a summary of our results. 

1. In the present context, phase separation involves the spontaneous formation of macroscopic domains of + and 
— as well as / and \ in the ULR Model || . This segregation is robust in that it survives at all temperatures 
T. Let us recall that most statistical systems which show phase separation at low T (or low noise-level, in 
nonequilibrium cases |) 10 1) lose this property at higher T or noise strengths. Certainly if one were to think in 



terms of energy and entropy, this would be the general expectation. Against this backdrop, a phase separation 
so robust as to persist at all finite T, and in a one-dimensional system at that, is quite unexpected. We suggest 
the name Strong Phase Separation (SPS) for this unusual phenomenon. 

The importance of SPS in the ULR model arises from the close relation of the latter to a physically realisable 
system of considerable current interest, namely, sedimenting colloidal crystals. Towards the end of this paper 
we suggest experiments which can be performed on fluidised beds of colloidal crystals to test some of the ideas 
presented in this work. 

2. The occurrence of SPS can be seen best in a certain limit in which the dynamics of the ULR model obeys the 
condition of detailed balance. In this limit, an energy function E can be constructed such that the steady-state 

probability of a configuration {e>"i,Tj_i} is proportional to exp — E({<7i, i }/T) . Although the dynamics is 
entirely local and involves rates of order unity, the emergent energy function E for the effective equilibrium theory 
involves interactions of unbounded range. As a result, E has a nonextensive (more properly, superextensive) 
character, which is how our model and those of |^,|| manage to get around the usual obstacles |ll[] to phase 
separation in one dimension. In our model, E has a simple interpretation: it is the energy of a collection of 
particles, viz. the ct^s, in a potential landscape built from {t^i}. The superextensivity is then a consequence 
of having potential energy wells whose depths scale with the system size. 
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Thermodynamic properties can be calculated in the strongly phase-separated state. In particular, the width of 
the interfacial region is found to vanish as T — > and diverge as T — > oo. 

3. Strong phase separation is a robust phenomenon, and persists even when the condition of detailed balance 
does not hold. This can be seen through arguments JEJ based purely on kinetics without recourse to an energy 
function: the transport of a + from one end of a + + + ...+ domain to a point a distance n away requires a time 
which grows exponentially with n, as n moves against the tilt field would be required. Thus a macroscopically 
phase-separated state would be expected to survive infinitely long in the infinite size limit. 

4. Although phase separation is inevitable in the ULR Model || and in the models of 0@]> the kinetics of domain 
growth is anomalously slow. The barriers that oppose the remixing of the macroscopically segregated state also 
inhibit the processes of diffusion that cause large domains to grow at the expense of smaller ones. These barriers, 
moreover, are produced by the dynamics of the model, not introduced ex machina in the form of quenched 
randomness. This results in intriguing ageing effects: for instance, the growth of domains is logarithmic in 
time, as has been verified in numerical studies in B. Further, in the detailed-balance limit, the decrease of 
E({ai, }) is logarithmic in time as well. Thus, despite the existence of a thermodynamic equilibrium state 
in the detailed-balance limit, a system which starts from a random initial condition has an extraordinarily 
difficult time reaching it. Such a system is best thought of as perpetually evolving, never in a truly steady state, 
sinking slowly into progressively deeper minima, in a manner which recalls the glassy state of the model of [ [L2[ . 

5. Arguments given in j^] already amounted to showing that SPS occurred in the ULR model. Specifically, it was 
shown there that the remixing of phase-separated domains in would always be opposed by barriers whose height 
diverged with the system size. The simulation results of j^], however, were complicated by the presence of a 
repulsion between adjacent + sites, which modelled interactions between charged colloidal particles. Increasing 
this repulsion beyond a threshold value led, in the numerical studies of Q, to an apparent loss ofphase separation. 
It is now clear, from the calculations reported in the present paper, that the observed remixing H] was a finite-size 
effect. 



D. Outline 

The rest of this paper is organised as follows. In the next section (|J) , we review the derivation || of continuum 
equations of motion for a crystalline array moving through a dissipative medium, and show how, at the linearised 
level, they lead to either a new class of "kinematic waves" [[l3| or an instability towards phase separation. Section 
H closes by presenting a simplified one-dimensional continuum model which retains all the essential features of the 
higher-dimensional problem. In section III, we use arguments similar to those connecting the noisy Burgers equation 
to the driven diffusive lattice gas to construct the LR lattice-gas model || whose long-wavelength limit has the 
relevant physics of the aforementioned one-dimensional continuum equations. We show, in a certain highly symmetric 
limit, that the Unstable LR model has a detailed balance property. In this limit we demonstrate Strong Phase 
Separation and calculate thermodynamic quantities. Further, we give arguments to show that SPS occurs in the 
entire parameter range of the ULR model. We argue that the coarsening of domains in the ULR model is ultraslow, 
with a characteristic length scale growing logarithmically in time. An analysis of a continuum model for SPS is the 
subject of section IV. Section V summarises our results and suggests experiments to test our predictions. 

II. CONTINUUM DYNAMICAL MODEL FOR A MOVING CRYSTAL 

A. Motivation 

The LR lattice-gas model || arose as a simplified description of the dynamics of a crystal moving steadily through 
a dissipative medium. It is therefore useful to review the construction of the continuum equations of motion for such 
a system. There are at least two physical situations where this dynamical problem arises: (i) the steadily sedimenting 
colloidal crystal mentioned above; (ii) a flux-point lattice moving through a thin slab of type II superconductor under 
the action of the Lorentz force due to an applied current. In (ii), the dissipation comes both from the normal core 
of the vortices and from disturbances in the order-parameter and electromagnetic fields in the region around the 
vortices. There is in principle an important difference between the sedimentation and moving flux-lattice problems: 
in the former, the disturbances produced by the moving crystal are carried to arbitrarily large length scales by the 
long-ranged hydrodynamic interaction, while in the latter, both electromagnetic and order-parameter disturbances 
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are screened and are thus limited to a finite range. A complete analysis of the sedimentation dynamics of a three- 
dimensional colloidal crystal thus requires the inclusion of the hydrodynamic velocity field as a dynamical variable. 
Instead, we consider an experimental geometry in which a thin slab of colloidal crystal (with interparticle spacing 
i 3> particle size) is confined to a container with dimensions L x , L z 3> L y ~ £ (gravity is along — z). The local 
hydrodynamics that leads to the configuration-dependent mobilities is left unaffected by this, but the long- 
ranged hydrodynamic interaction is screened in the xz plane on scales 3> L y by the no-slip boundary condition at the 
walls. The model equations (0) in dimension d = 2 apply to such a system. 



B. Constructing the equations 

Our construction of the equations of motion ignores inertial terms, which is justified both for the confined colloidal 
crystal and, except at very low temperatures [jl4|| , for the flux lattice. Rather than keeping track of individual particles, 
we work on scales 3> £, treating the colloidal crystal or flux lattice as a permeable elastic continuum whose distortions 
at point r and time t are described by the (Eulerian) displacement field u(r,i). In general, the equation of motion in 
the completely overdamped limit has the form velocity = mobility x force, i.e., 

^u = /i(Vu)(KVVu + F + f). (2) 

In (^), the first term in parentheses on the right-hand side represents elastic forces, governed by the elastic tensor K, 
the second (F) is the applied force (gravity for the colloidal crystal and the Lorentz force for the flux lattice), and f is 
a noise source of thermal and/or hydrodynamic origin. Note that in the absence of the driving force F the linearised 
dynamics of the displacement field in this overdamped system is purely diffusive: dtu ~ V 2 it, with the scale of the 
diffusivities set by the product of a mobility and an elastic constant. All the important and novel physics in these 
equations, when the driving force is nonzero, lies in (j,, the local mobility tensor, which wc have allowed to depend on 
gradients of the local displacement field. The reason for this is as follows: The damping in the physical situations we 
have mentioned above arises from the interaction of the moving particles with the medium. A dynamical friction of 
this kind will in general depend on the local arrangement of particles |j 
the symmetry of the mobility tensor will thus reflect the symmetry of the underlying lattice. If the structure in a 
given region is distorted relative to the perfect lattice, the local mobility will depart from its ideal structure as well. 
Deviations of the structure from the perfect crystal are described by the full distortion tensor Vu |l(| rather than 
its symmetric part, the strain, since we are not in a rotation-invariant situation. We further make the reasonable 
assumption that the mobility can be expanded in a power series in the distortion: 

Ax(Vu) = Mo + AVu + C>((Vu) 2 ), (3) 

where fio is the mean macroscopic mobility of the undistorted crystal. 

For a d-dimensional crystal driven steadily along the z direction, assuming isotropy in the d — 1-dimensional "_L" 
subspace normal to z, but not under z — * —z, (||) and (||) lead directly to 

u'i. = Ai<9 2 u ± + \ 2 V±_u z 

+C(VVu) + 0(VwVu) + f x , (4a) 
u z = A 3 Vi.ui + Xid z u z 

+C(VVu) + OiVuVu) + f z , (4b) 



15 . Even for a perfect, undistorted lattice, 



where the constant drift along z has been removed by shifting to the mean rest frame of the crystal. The terms that 
are manifestly most important at small wavenumbers, at least within a linear description, are the linear, first-order 
space derivative terms. These terms arise from (Q) and (|J) via the leading distortion-dependence of the mobility 
tensor, multiplied by the driving force F. The coefficients A^ [as well as those of the 0(VuVu) terms, as can be 
seen from (||) and (Q)] are thus proportional to F, and the corresponding terms are therefore present only in the 

driven state. At small enough wavenumbers ( ^ F/K where F is the magnitude of the driving force density and K 
a typical elastic constant), these terms dominate the diffusive terms coming from the elasticity. The terms of this 
type in (|4a|) tell us that a tilt (a z derivative of a _L displacement or a 1 derivative of a z displacement) leads to 
a lateral drift, and those in ( pb| ) imply that the vertical settling speed depends on the compression (or dilation). 
Since the system is not invariant under rotations, there are no grounds for insisting that Ai = A2 or A3 = A4. f 
is a spatiotemporally white noise source containing the effects of thermal fluctuations as well as chaotic motion due 
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to the hydro-dynamic interaction |17],|l8| . The reader will note that the form of the diffusive second derivative terms 
and the distortion-dependence of the mobility beyond linear order has been left rather general. This is because even 
for d = 2, as can be seen by exhaustive listing, symmetry under permits, all told, in (|4a|) and 

(^h|), ten terms (this counting was wrong in j^j) bilinear in Vu and six linear second derivative terms, with as many 
independent coefficients. It is clearly difficult to make very useful general statements about a problem with so many 
phcnomenological parameters so we restrict ourselves, in the next subsection, to a linearised description to lowest 
order in gradients. We will return to the effects of nonlincaritics in later subsections. 



C. Mode structure 



If we retain only terms linear in the fields and work only to leading order in wavenumber, then the relation between 
frequency ui and wavevector k implied by (|j) is 



— I 



(Ax + \ A )k z ± ^(Ai -A 4 ) 2 fc2 + 4A 2 A 3 fc2 



(5) 



The dispersion relation (^J) has a wavelike character in all directions if A2A3 > 0. For A2A3 < 0, while it is still wavelike 
for k± = 0, it has a growing mode to oc — ik for k z <C k±. 

Linearly stable case — kinematic waves: The wavelike modes are the generalisation, to the case of a moving lattice, 
of the kinematic waves which Lighthill and Whitham [jl3| discussed in the context of traffic flow and flood movements. 
The important difference in the present case is that the waves propagate not only along but also transverse to the 
direction of drift. Some remarks towards a more complete consideration of their dispersion relation, including the 
effects of nonlinearities, may be found in the context of a one-dimensional reduced model in j5|. 

Linearly unstable case — clumping: In the case A2A3 < 0, for wavevectors pointing outside a cone around the z 
axis, the system is linearly unstable, as already noted in [^J: small perturbations grow, with a growth-rate which is 
linear in their wavenumber. Whereas the linearised treatment cannot give detailed information about the final state 
of the system, we expect the growing mode to appear as a clumping and tilting of the colloidal crystal, with material 
concentrated at the bottoms of the tilted regions. The wavevector of the inhomogeneity will be mainly normal to the 
scdimenting direction. 

The remainder of this paper is directed towards a more detailed understanding of the statistical mechanics and 
dynamics of macroscopic clumping. Our studies are based mainl y on t he simplified one-dimensional lattice model of 
II . The construction of the lattice model is reviewed in section HI A: its origins lie in a reduced, one-dimensional 
version of equations (0) which we now present. 



D. A one-dimensional effective model 



We saw above that the equations of motion for a moving lattice contained terms of a qualitatively new form, not 
present in the equations of a lattice at equilibrium. To linear order, these were the {Xi} terms in (^), which are 
proportional to the driving force, and of lower order in gradients than those arising from the elasticity of the crystal. 
The effects of the linear instability for A2A3 < thus cannot be mitigated by including the diffusive terms arising from 
the linear elasticity. To see what final state, if any, emerges from the initial unstable growth in the case A2A3 < 0, 
we must go beyond a linear treatment. Even in the stable case A2A3 > 0, the combined effects of nonlinearities and 
noise could result in effective dispersion relations for long wavelength modes which differ qualitatively in their form 
from those predicted by the linear theory. However, including nonlinearities, diffusion and noise, as we remarked in 
the previous subsection, introduces an enormous number of phenomenological parameters. We note instead that the 
important new physics of (Q), namely, the wavelike (stable case) or growing modes (unstable case), arises from the 
coupling of the vertical and horizontal displacement fields, for excitations with wavevector transverse to the direction 
of mean drift, while the modes with wavevector along z play a relatively minor role. This suggests that much can be 
learnt from a model in one space dimension, the x direction, corresponding to the _L direction of (^) , but retaining a 
two-component displacement field u — (u x ,u z ). The symmetry X ^ X ^ t-^x ^ "^x 

then yields, to bilinear order in 

fields and leading orders in gradients, the equations of motion 

u x = X 2 d x ii z + r )\d x u x d x u z + Did x 2 u x + f x (6a) 
u z = X 3 d x u x + j 2 (d x u x ) 2 

+l 3 (d x u z ) 2 + D 2 d x 2 u z + f z , (6b) 
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which have, in addition to the {A^}, three nonlinear coupling parameters {7^} (also proportional to the driving force F), 
two diffusivities {A}, and gaussian spatiotemporally white noise sources fi, i = x,z, with zero mean, and variances 
N x , N z : 

(fi(0,0)fj(x,t)) = 2N i 5 ij S(x)5(t). (7) 

If {ji}, {Di} and {/,-} are set to zero, we recover the continuum limit of the equation derived by Crowley || for the 
dynamics of the small transverse and longitudinal displacements of a collection of hard spheres of radius a, prepared 
initially in a horizontal, one-dimensional periodic array with spacing d, settling vertically in a highly viscous fluid, with 
the hydrodynamic interaction cut off at the nearest neighbour scale. The correspondence is A2 = —A3 = — (3/4)a/d, 
in units of the Stokes settling speed of an isolated sphere. Crowley's calculation can be extended beyond linear order 
to give {7^}, but the elastic forces and the thermal fluctuations that give the DiS and fiS are absent in his model. 
The diffusion and nonlinear terms in (||) are identical in structure to those in the Erta§-Kardar (EK) models for the 
fluctuations of drifting lines jl9|,^0|], with u x , u z replaced by their variables h±, hu in jl9| or R±, Rn in ^o|. The 
EK models, however, as a result of a larger symmetry (independently under (i) x — » — x and (ii) R± — » —R±_ or 
h± — > — hx) lack the linear first spatial derivative terms (the Ai terms) of (||). Such linear terms can however be 
induced through the nonlinear terms, in [ p~9|j20| by constraining the ends of the line (polymer) to be at fixed mean 
separation normal to the drift direction, so that ) ^ 0. The related coupled- interface model of Barabasi |2l]] 
has ann — x symmetry and thus also lacks the Ai terms of (||). These models are thus not relevant to the case of 
greatest interest to us here, namely, the unstable case A2A3 < of (^|). 

In the unstable case, within a linear treatment, the concentration d x u x and the tilt d x u z grow without bound p2j. 
Physically, since real colloidal crystals are made of impenetrable particles, and since the elasticity of the lattice will not 
tolerate arbitrarily large shear-strains, the description implicit in (^|) of small distortions about a perfect lattice must 
break down in conditions of unstable growth. It is best, therefore, to work from the outset with naturally bounded 
variables for the concentration and tilt. To this end, we first pass to a description in terms of the concentration 
fluctuation field 

(8) 



and the tilt field 



r{x,t) = -^. (9) 

Then (^]) can be rewritten in the "conservation-law" form 

cr = A 2 5 x r + jid x ((TT) + Did x 2 (T + d x f x (10a) 
f = X 3 d x a + j 2 d x (ct 2 ) 

+73^(r 2 ) + D 2 8 x 2 t + d x f z . (10b) 

As stated above, a and r should be bounded; what matters on large length scales is only whether the local concentration 
is large or small compared to the mean, and whether the local tilt is "up" or "down" . Accordingly, we construct a 
description in the next section in which the concentration and tilt fields of ( [To|) are replaced by Ising variables evolving 
under a spin-exchange dynamics designed to mimic the most important aspects of (|l0|). A continuum model which 



incorporates saturation is presented in section IV 



III. STRONG PHASE SEPARATION IN A LATTICE MODEL 

In this section, we introduce the notion of Strong Phase Separation in connection with the LR lattice model, which 
describes two coupled species of spins on a lattice, with simple evolution rules which mimic the coupled dynamics 
of the density and tilt fields. This coupled-spin problem is too difficult to solve for the dynamics or, indeed, for 
the steady state for arbitrary values of parameters. However, for th e sym metric case of half filling of both species, 



and a special relation between coupling constants, we show (Section IIIB) that the condition of detailed balance is 



satisfied with respect to a Hamiltonian Ti. with long-ranged interactions. In turn, this allows for a characterization 



of the steady state of the system. In Section [II C , we show that at zero temperature T, the system exhibits phase 
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separation. Moreover, we calculate thermodynamic properties and show that the phase separation survives at all 
finite temperatures, which is why we call this phenomenon Strong Phase Separation (SPS). The occurrence of SPS is 
linked to the long (actually infinite) range of the interactions in TC, which results in the energy being superextensive 
(proportional to L 2 rather than L). We emphasise that this happens although the underlying dynamical model is 
entirely local, with finite, bounded rates. In section HID, we show that this unusually robust phase separation sets 
in anomalously slowly, with domain sizes growing as the logarithm of tim e. Th e survival of SPS away from the 
detailed-balance limit is discussed through a kinetic interpretation in section III E| . 



A. The LR Lattice Model 



^From the study of driven diffusive systems, it is well known that hydrodynamic behaviour can be recovered from 
the large-distance long-time behaviour of simple lattice gas models evolving by stochastic dynamics [^3| . An example 
of such a model is the asymmetric exclusion process, in which particles on a lattice perform biased random walks 
subject to the constraint of no more than one particle per site; in the limit of large separations and time, density 
fluctuations are described by the Burgers equation with an additional noise term. An advantage of a lattice gas 
description is that nonlinearities are incorporated implicitly in the nature of the variable - for instance, a (0, l)-valued 
occupation variable incorporates the effects of exclusion. 

Are there simple lattice gas models which capture the essential features of coupled density-tilt dynamics of the type 
discussed in the previous section? Any such lattice model must, of course, involve two sets of variables - say {<Ji} 
and {r,} - which are discrete versions of density and tilt fields and which evolve by rules which mimic the physics of 
sedimenting lattices. There are two crucial features of the a — r dynamics of Eq. ([To]): first, that both a and r fields 
are conserved so that their time derivatives involve the divergences of currents; and second that the local field which 
guides the cx-current has a term which is proportional to r, and vice versa. Accordingly, we define || a lattice model 
which incorporates just these effects. Consider a one-dimensional lattice made of two interpenetrating sublattices S 
(i = 1, 2, 3, • • • , N) and T (i = |, |, • • • , N — h). Place Ising variables (Ji = ±1 at every site of S, and r i+ i = ±1 on 
every site of T. We take Oi — 1 if there is a particle at site i, and <7; = —1 if there is no particle, while Tj = 1 or 
— 1 denotes the two possible values of the local tilt. The dynamics involves exchange of adjacent spins <Ji and c^+i at 
a rate which depends on the intervening spin while the rate of r-spin exchanges depends on the intervening a 

spin, i.e. we have Kawasaki |24|] dynamics, with hopping rates which depend on the local value of the other species. 
The probability P(C) that the system is in a configuration C = {{<?i}, { T i-i }) evolves through the master equation 

^jp- = W(C n , n+1 ^C)P{C n , n+ i)-W(C^C n , n+ i)P(C). (11) 

(n,n+l) 

Here (n, n + 1) on the right hand side (with n — i, 1, ~, 2, • • •) labels transitions which involve pairwise interchanges 
of neighbouring er's (cr^ «-> <Ji+\) and r's (t^i «-» and configuration C„ in +i differs from C only through the 

interchange of spins on site n and n + 1. The most general such model would involve the 8 distinct transition rates 
listed in Eq. (1). For a left-right symmetric system, we have D — D'; a — a'; sgn(b) = sgn(b'): this defines the LR 
model ||. 

In the interest of defining a minimal version of the LR model, we also impose the further restrictions E = E' , 
b = b' . The rates of the minimal model may be written compactly as 

ar i+ i 

W{a l <-> cr i+1 ;r i+ i) = D jT^fa ~ 



W(T i _ i ~T i+i ;< H )=E+^(r i _ i -T i+i ). (12) 

The evolution rules can be stated as follows: If a is positive, a particle tends to move downhill, and a hole uphill. If 
b is positive, a local peak (A) tends to transform into a valley (V) if a particle resides on it, while local valleys tend 
to become peaks in the presence of holes. Changing the signs of a and b reverses these tendencies. As a result, the 
nature of the steady state is sensitive to the sign of a = ab. As we will see below, if a is positive, the exchanges of 
a and r spins in Eq. (|l2] ) act in concert to promote segregation of both species of spins, ultimately resulting in a 
phase-separated state. This is the unstable case of the LR model - the case of primary interest in this paper. By 
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contrast if a is negative, 'easy' a and r moves produce opposing tendencies, and hence result in a fluctuating but 
on-average spatially homogeneous state - the stable case of the LR model. The calculations of Crowley (f| for settling 
arrays of hydrodynamically interacting spheres and the discussion in || make it clear that for sedimenting colloidal 
crystals it is the "unstable" case that applies. 

The other important parameters in the model are the magnetizations M a = £i<7,;/iV, M T — SjT i+ i/iV, both of 
which are conserved by the dynamics. 



B. Symmetric Case: Hamiltonian and Detailed Balance 



We now consider the symmetric case of the LR model, which is defined by the vanishing of the magnetizations 

M c = M T = 0, (13) 

and the following relationship between coupling constants in ( |l2|) 



b a 
E ~ D' 



(14) 



Since E, D > 0, it is clear that (14) is a special case of the unstable LR model. We show below that when conditions 
( |T3| ) and ( p^ ) are met, it is possible to find a Hamiltonian TL such that the condition of detailed balance is satisfied 
with invariant measure exp(— PTC). 

Since the motion of a particles is determined by the local tilt r, we may think of the a particles as moving in a 
potential landscape provided by the r's (Fig. With this in mind, we define the height at site k by 



(15) 



With periodic boundary conditions (ctn+i — Ci] T N+i _i = r i _i), the zero-net-tilt condition M T = implies hjsr+k 
hk- We associate a potential energy proportional to hkCk with site k, and write the Hamiltonian 



N 

H = h k {r}(7 k 

fe=i 



(16) 



to describe the total energy of the a particles in the landscape derived from the r particles. 




□ODODODODOD9D* ■•■•■•■OBOBOBO BO BO 

FIG. 1. The phase-separated state of the 1-d lattice model at T=0 is shown, a and r variables are shown as circles and 

squares respectively, with a, r = +1( — 1) shown filled (empty). The configuration of the corresponding height model is also 
shown. Interfaces between a — +1 and —1 are located at A and A' and those between r = +1 and —1 are at B and B' . 
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In view of the symmetric role played by er's and r's in the Symmetric Model, we may equally ask for the potential 
energy of t particles in the landscape provided by the a particles. The corresponding Hamiltonian is then 

N 

Q = ^st-^W^-i ( 17 ) 
k=i 

where the height g is given by 

fc 

j*+iW=E ff i- ( 18 ) 

When the heights hk and 9k-\ are written out in terms of r^'s and er^'s respectively, the Hamiltonians TL and Q are 
seen to involve very nonlocal couplings: 

N k 

fc=l J=l 



N k-1 

e = e EE^ r H- (20) 

We observe that 

W + G = eM a M T (21) 

and, since each of M a and M T vanishes in the symmetric case owing to the zero-tilt condition, we have Ti = —Q, 
Thus the Hamiltonians corresponding to the two pictures i.e. a particles in a r-landscape or vice versa are completely 
equivalent. We will mostly use Ti for further work. 

We now show that the steady state of the Symmetric model defined by Eqs. ( |l3| ) and (|l4|) satisfies the condition 
of detailed balance and that the stationary measure is given by e _/3W where (3 is the inverse temperature T -1 , with 
/3e given by below. To this end, let us ask for the changes in energy AE(ai «-» <Ji+i) of Ti when spins cr, and 
are interchanged, and A£ , (r i _ 1 / 2 Ti+1/2) m W when spins and T i+1 / 2 are interchanged. For i ^ JV, it is 
straightforward to see that 

A£(cTi <-> o-j+i) = eT 4+ i(cri - cr i+ i) (22) 



AS( Ti _i ~T i+i ) = ea i (r i+i -T i _ i ) (23) 

In fact, Eqs. (pffi) an( i (H) are valid for i = iV as well, as can be verified on recalling that <jn+i — t n+± = r i/2 
and using the zero-tilt conditions M a = M T — while computing energy changes. 

Consider the configuration C CTil0 - j+1 obtained from a configuration C on exchanging two neighbouring a spins - an 
elementary kinetic move in the model. The condition of detailed balance is then 

W{C^C^ +1 ) = Mss(C CTt ^ +1 ) 

W{C ai , ai+1 - C) mss(C) 1 j 
where /is's(C) is the steady-state measure for configuration C. To verify that 

/lss (C)=e-^, (25) 

we use Eqs. ( |l2| ) and (^) to obtain 

SiS-'** (26) 

where we have defined X, = (<Zj — cr; +1 ). Noting that Xj = ±1, we see that Eq. (|2^) is satisfied provided 







0e 



2 \D — a 



(27) 



In order for the measure to be valid under interchanges of adjacent r's iji-i T i+ i), similar reasoning leads to the 
condition 



(28) 



and ( P8[) are consistent. Thus the condition 



In the Symmetric case of the LR Model, Eq. ( |14|) holds, and so Eqs. 
of detailed balance holds with the equilibrium measure (25). 

It is appropriate to recall that the three-species model of Evans et al. ffjj also obeys the condition of detailed balance 
in the symmetric case. There too the Hamiltonian has infinite ranged interactions, but does not have as transparent 
an interpretation as fll6|). 



C. Symmetric Case: Thermodynamic Properties and Strong Phase Separation 



Since the condition of detailed balance holds in the symmetric case of the minimal LR model, the steady state 
corresponds to the thermal equilibrium state with Hamiltonian Ti. The thermodynamic properties of the system can 
be found, in principle, using equilibrium statistical mechanics. A calculation can be carried out in the grand canonical 
ensemble in the limit A — » oo. The resulting state exhibits Strong Phase Separation. 

The Hamiltonian H (Eq. ©) describes spins <Jk in a site-dependent magnetic field ehk, which is itself a dynamical 
variable. Equivalently, in the lattice gas description (associating an occupation variable nk = |(1 + crk)), it describes 
particles with a hard core constraint in a potential well of depth ehk- The ground state of H is obtained by arranging 
the r spins (which determine the heights hk) so as to form as deep a potential well as possible, and then arranging 
the c-particles at the bottom of the well (Fig. [if). A spin configuration which corresponds to this choice is 

T fe _ 1/2 = -1 for k = !,-■■, N/2 

A 

= 1 for k = — + 1,---,N 

a k = 1 for k = A/4, • • • , 3A/4 

= -1 for k = 1, • • • , A/4 - 1 and k = 3A/4 + 1, • • • , A. (29) 

Each spin species exhibits complete phase separation in this ground state. The ground state energy is straightforward 
to compute, and we find 

eA 2 

Eg * - — ■ (30) 

Notice the quadratic dependence of Eq on A, which is an outcome of the infinite-ranged interactions in Ti (Eq. [30| ). 
As explained below, this unusual superextensive behaviour of the energy is ultimately the feature responsible for the 
phenomenon of Strong Phase Separation, namely the continued existence and stability of the phase separated state 
at all finite temperatures. 

At T = 0, phase separation is complete and there is a sharp boundary between regions of positive and negative 
spins of each species. Let A and A 1 be the locations of the T = interface between regions with a = 1 and a = — 1, 
and let B and B' be the locations of interfaces separating regions with r = 1 and r = — 1 (Fig. |l|). The effect of 
raising the temperature to a finite value T is to smear out the interfacial zones around A, A' , B and B' (Fig. ||). To 
address this quantitatively, let us turn to the evaluation of thermodynamic properties. 
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FIG. 2. Typical configurations of the Unstable LR model at T / are shown. The meaning of the symbols is as in Fig. 

Phase separation persists, but there are particle-hole excitations of both species near the corresponding interfaces. 

The calculation can be carried out most easily in a grand ensemble in which the total magnetizations M a and M T 
are not held fixed. The corresponding grand partition function is 

Z= e-K n - E <* 

= e^S*^^-^) (31) 

MM 

where <r° denotes the value of <Jk in the ground state. The key observation that allows the calculation to be performed 
is that near the u interfaces A and A', the field hk({r}) is essentially fixed at its T = value h? k ; deviations are of 
order exp[— (3eN/A\ as explained below, and so are utterly negligible in the thermodynamic limit. Likewise, in the 
vicinity of the r-interfaces B and B', the a spins are frozen to their T — values, and so <7fc({c}) = g®- To proceed, 
let us divide the system into four equal parts Ra, Rb, Ra', Rb', where region Ra consists of the N/A spins of each 
of the two species centered around A. Other regions are defined similarly, centered around B, A' and B' . Evidently, 
with negligible error we may set hk({r}) = hi in regions Ra and Ra 1 , and set 3fe({c}) = g® in regions Rb and Rb' ■ 
The partition function Z can then be written as the product of 4 terms Za, Zb, Za' , Zb> , where for instance 

Z A = y_. e ^ keR A kK k ' (32) 

z B = j2 e ~ 0eEk ' RB9l+i{Tk+i ~ T ° k+i) - ( 33 ) 

M 

Each of these factorizes into single-site partition functions, and can be evaluated straightforwardly. Recalling that h® 
varies linearly with k near the T = interface location Ua, we find 

Za = II (l + e~ 2Mk ~ kAl )- (34) 

k£R A 



In the thermodynamic limit, we obtain 



Z A = K(e~ 2 ^) (35) 



where K{y) = rifcL-oo(l + V 1S a generating function that arises in the theory of partitions p5| . Evidently, each 
of Zb, Za> and Zb> equals the same quantity as well, so that Z — [K(ex-p(— 2/3e))] 4 . 

It is worth pausing to comment on the unusual size dependences of various quantities. The ground state energy 
Eg is proportional to -/V 2 , a superextensive dependence. This has its origin in the infinite-ranged interactions in 7i. 
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Further, with energies measured from the ground state value, the partition function approaches an TV-independent 
limit. Thus the total change in free energy and entropy away from T = remain finite in the thermodynamic limit 
i.e. they are not extensive. This reflects the fact that the only effect of raising the temperature is to broaden the 
interfacial region between phases, which essentially affects only a finite number of sites. 

In fact, an explicit calculation of the broadened interfacial profile can be carried out in the grand ensemble. For 
instance, near A we have 

(cr fe ) =twh^ehl (36) 

where hk = (k — k A )- We see that (<jf.) deviates substantially from 1 only in a region where \/3ehk\ ^ 1, or 

\k-k A \<T/e. (37) 

For sites k such that \k — k A \ 3> T/e, the deviation from ±1 is ~ 2exp(— 2(3e\k — k A \) which vanishes rapidly. We see 
that the primary effect of temperature is to smear out the interfaces. The formation of 'droplets' far from the interfaces 
is prohibitively costly in energy, and hence the probability dies down exponentially. Recalling that the separation of 
the two a = 1 — ► — 1 interfaces is N/2, in the thermodynamic limit N — > oo, we see that only a vanishing fraction of 
spins (those close to the interfaces) deviate from values arbitrarily close to 1 and —1. In this sense, phase separation 
remains complete and cannot be effaced at any finite temperature T, i.e. we have Strong Phase Separation. 

These results obtained in the grand ensemble provide a qualitative, if not quantitative, guide to the thermodynamic 
properties of the system in which M a and M T are held fixed. The customary equivalence between ensembles is not 
obviously valid any longer, as particle-hole excitations are essentially confined to a finite region of width proportional 
to T, which does not increase as N — > oo. Thus the difference between observables calculated in the two ensembles 
is expected to remain of order unity, and not die out in the N — > oo limit Interestingly, the calculation of the 
partition function, though not the profile, has been carried out for the three-species model within a constant-species- 
number ensemble 0. 

The stability of the strongly phase-separated state can also be understood in terms of kinetics. In the ground state 
arrangement of Fig. |], each a spin finds itself in a uniform field produced by the r spins. Consider moving a spin over 
a macroscopic distance - say a a = +1 spin from A 1 to A, via B. The movement from A' to B may be viewed as an 
activation process as the spin in question has to overcome a potential barrier of magnitude eN/4 to reach B; beyond 
that, in the region BA, the motion is ballistic as the r-induced field helps it along. The rate-limiting step is thus the 
A' — > B activation. At temperature T, the relevant time scale is of the order oitcB ~ exp(e./V/4T) which diverges 
rapidly as N — > oo. Thus, in the thermodynamic limit, a rearrangement of the SPS state is not possible; the only 
effect of the temperature-assisted motion is to move a few a = 1 spins near the interface into the er = — 1 rich region 
and vice versa, but such penetration does not proceed far in view of the restoring fields. Defining the penetration 
depth Afc as that over which the activation time falls by a factor of 1/e, we estimate Ak = T/e, in agreement with 
Eq. ( |37| ) which was based on the spatial decay of the interfacial profile. 



D. Coarsening 

Now imagine that Fig. [l] represented half the system, and that the other half was identical in structure. This 
would amount to a system that had phase-separated into four macroscopic domains, each of size N/8. For this state 
to proceed towards full phase separation, the two + domains, each at the bottom of a valley, must merge. The 
rate-limiting step can again be taken to be the movement of a + from the edge of an all — h region to the top of a hill, 
i.e., a distance N/8. Once this comes to pass, the two domains of length N/4 will rapidly merge to give one domain 
of length N/2. The time for t his, w hich is the time for complete phase separation for a system of size N, can be seen 
from the argument in section [II C| to scale as exp(eiV/8T). This tells us that the characteristic domain size grows 



logarithmically in time, as stated in section [|. 

The time required for the reverse process (from a 2-domain to a 4-domain state) scales as ~ exp[(iV/4T)], which 
is overwhelmingly larger than the 4 — > 2 coarsening time. This is true at all scales, and the transition from a 2n- 
domain state to one with n domains is much more rapid than the reverse. Thus the transition from a statistically 
homogeneously mixed state to the equilibrium phase-separated state is irreversible, even though it occurs slowly. 

The coarsening process was studied [Q both numerically in the 3-species model and within a mean-field approxi- 
mation for a related 'toy' model. The typical domain size was found to grow logarithmically in time. The arguments 
given above are consistent with this. 
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E. Non-symmetric Case 



We now address the nature of the steady state for arbitrary values of M a and M T . Away from M a = M T = the 
problem is no longer described in terms of the equilibrium state of a long-ranged Hamiltonian; nevertheless we will 
argue below that the system continues to exhibit Strong Phase Separation. 

It is useful to define x and y as the density of up spins of the a and r types. We have x = |(1 + M a /N) and 
y = |(1 + M T /N). If x and y are small enough that 2x + y < 1, the steady state is of the type shown in Fig. |]a, with 
each of the a and r species showing phase separation, but with basically no spatial overlap of the a = 1 and r = 1 
regions. A useful way to characterize this state is through the sequence of interfaces, viz. A ■ ■ ■ B' ■ ■ ■ BA' ■ ■ •, where 
• • • denotes a macroscopic stretch of the system. Here A{B) separates an up-spin region of <j(t) spins on the right, 
from the corresponding down-spin regions, while A' and B' separate the opposite regions. Trial states of the type 
A ■ ■ ■ B' ■ ■ ■ B ■ ■ ■ A' ■ ■ ■ are seen to approach the non-overlapping state on a time scale of order t* where In t* is less 
than but of the order of the smaller of eNx/T and eNy/T. Once the nonoverlapping steady state has been reached, 
a and r spins can still be cycled around by activation processes across A' A and A' A respectively (Fig. |^), but such 
cycling around does not change the character of the state. 




(a) (b) 

FIG. 3. Typical configurations are depicted away from the symmetric case in the limit of vanishing noise, (a) If the fraction 
of a = +1 and r = +1 spins is low enough, interfaces A' and B coincide, (b) If the fraction of a = +1 spins is high enough, 
the interface B' lies halfway between A and A' . 

Now consider increasing y, keeping x fixed. The number of spins in the stretch between B 1 and A is 7V(1 — x — y), 
and once this drops below Nx, the predominant activation process occurs over this stretch. Thus the no-overlap state 
of Fig. H^, is unstable towards a state of the type shown in Fig. ||b, once 2x + y exceeds unity. In this steady state, 
activation processes in a finite system lead to small currents of a and r spins, of magnitude 

J a = ai cxp(- el ab 1 /T) - a 2 cxp(-ef j4 'B'/ r ) ( 38 ) 



J T = a 3 cxp(-e£ B A/T) - a 4 exp(~e£ B 'A/T) (39) 

where 01,02,(13,04 are prefactors of order unity and Iab 1 is the separation of interfaces A and B' , and other Ps 
are defined similarly. Since the difference £a> b' ~ £ab' is positive and grows proportionally to N, we may drop 
the second term on the right of Equation (|38|). In steady state we must have J CT = J T , which then leads to (a% + 
a 4 ) cxp(~ e£ AB > /T) ~ a 3 exp(-e£ B A/T) or 

£ab' = £ba + terms of order unity. (40) 

Thus, A is very close to the halfway position in stretch BB' . The overlapping stretch BA' is a fraction S = ^ (2x+y— 1) 
of the whole length. On setting x — y = h , we recover 5 = 1/4 in agreement with the results of the equilibrium 
analysis of Section C. 

Analogously, keeping y fixed and increasing x we conclude that for 2y + x > 1, the steady state is 
A - ■ ■ B ■ ■ ■ A' ■ ■ ■ B' ■ ■ ■ with B' between A' and A and an overlapping stretch of length ^N(2y + x — 1). Finally, 
under x — > (1 — x), y — > (1 — y) we arrive at the condition for overlap of negative spins. 
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In short, strong phase separation persists even away from the symmetric point of the LR model. In general, two 
types of steady states, both phase separated, are possible as depicted in Fig. ||. In the overlapping case, there is 
generally a current in a finite system, but this vanishes exponentially with system size. While we have explored the 
effects of deviating from the symmetric case by moving away from the half- filling condition Eq. (13), without altering 
the condition (14) on the rates, another way to make the system nonsymmetric is to violate the latter condition. We 
have not explored this in detail but expect that the phenomenon of SPS will persist in this case too so long as ab > 0. 



IV. DETAILED BALANCE AND STRONG PHASE SEPARATION IN THE CONTINUUM LIMIT 

The continuum model of section in the case A2A3 < in equation (^J), is linearly unstable. One way to deal with 
this instability is to resort, as we have done above, to a lattice model in which the variables are naturally bounded. 
An alternative way is to ask what nonlinear terms added to @ for A2A3 < would arrest the unstable growth [p2[ . 
To do this, we work in the detailed-balance limit of the lattice model, start with the Hamiltonian ([l6|), and construct 
the corresponding continuum Ginzburg-Landau free-energy functional. We shall see below that this functional will 
give rise to dynamical equations with the same linear instability as in (^|) with A2A3 < 0, but containing nonlinearities 
which prevent unbounded growth. 

The derivation is straightforward, as the condition of detailed balance allows us to proceed as in any equilibrium 
statistical mechanics problem. The Ginzburg-Landau free-energy functional F[o~, t] for our syste m, i.e., the effective 



Hamiltonian for a description in terms of the coarse-grained fields {<r(x), t(x)} of section II D , may be written as 
U — TS where U is the energy (Jl^) in the continuum limit, T the temperature, and S the entropy obtained by summing 
over all microscopic configurations {eri,Ti} subject to a fixed coarse-grained configuration {a(x), t(x)}. Since o~i and 
t, are Ising variables, S can be found from a standard Bragg- Williams construction. Thus, 

F[u,t]=€ [ dx [ Ax'(j{x)t{x') 
Jo Jo 

<m 1 -L nn 1 — rn 1 — ry) 

(41) 



m—T (x) : <j{x) 



(i 



1 + m , 1 + m 1 — m , 1 — m 
In 1 In 



where x is measured in units of the lattice spacing and is hence dimensionless. The partition function is then 
f[da}[dr]exp(-F/T). 

Then the usual, purely dissipative, conserving time-dependent Ginzburg Landau equations of motion generated by 



i.e., 



SF 

d t a = A r7 d 2 x —+r, a (42) 



and likewise for r, turn out to be precisely 

d t <j = A CT (T^tanh" 1 cr + ed x r) + n a ; 

d t r = A T (T^tanh" 1 T - ed x a) + rj r . (43) 

Here A CT , A r are mobilities and rj a , rj T are noise sources with variances proportional to the corresponding mobilities. 
It is evident that equations (^3|) and ( |l0| ) are identical in the linearised limit, if we make the identification A2 = 
A CT e, A3 = — A r e. This corresponds to the linearly unstable limit of (^), in consonance with fact that the detailed 
balance limit of the lattice model was derived in precisely that case. 

We should thus be able to gain some insight into SPS by looking at the steady states of ([43]). The simplest of these 
are the zero current states, which satisfy 

d x P - tanhQ = 0, 

dzQ + tanhP^O, (44) 

where 

Q ee tanh -1 <?!>, P = tanh" V (45) 

The spatial development of P and Q with respect to x is like a Hamiltonian dynamics, conserving the "energy" 

E(P, Q) = ln(cosh P cosh Q) (46) 

This leads to closed orbits in the P — Q or ip — <j> plane, i.e., regions of large ip and small <j> followed by the opposite. 
These are spatially multidomain states which will not evolve further in the absence of noise. 
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V. SUMMARY AND DISCUSSION 



A. Summary 

In summary, we have constructed continuum and lattice models to describe the physics of steadily sedimenting 
colloidal crystals or, more generally, of a crystal driven through a dissipative medium. The models display two 
broadly distinct types of behaviour, termed "stable" and "unstable", depending on the sign of a parameter. We 
have concentrated on the unstable case and shown, through a mapping to a one-dimensional lattice model, that it 
always displays phase separation, a phenomenon which we call strong phase separation. This phase separation and 
the fact that it persists at all temperatures can be understood, in general, in terms of barriers to remixing which 
grow with system size. The barriers are erected by the system in the course of its evolution, and result in domain 
sizes growing as the logarithm of the time. In a particular limit, the detailed balance condition holds, allowing us to 
write the steady state distribution in the equilibrium form exp(— [3TL), and to calculate density profiles exactly. Here 
TL involves long ranged interactions even though the model has strictly local dynamics. This long-ranged character of 
interactions in 7i is responsible for the phase separation in this one-dimensional system, and the fact that it persists 
at all temperatures. 

B. Experimental tests 

Finally, let us turn to the possibility of testing our results in experiments. We have demonstrated Strong Phase 
Separation in a one- dimensional model system. It seems highly likely, therefore, that the same phenomenon will take 
place in the experimental systems which inspired our model, namely, steadily sedimenting crystalline suspensions in, 
for example, the two-dimensional geometry described in section |. A good candidate system is a charge-stabilised 
crystalline array of polystyrene spheres with radius in the micron range. The lattice spacing of the crystal should 
be neither so large that hydrodynamic effects (proportional to the ratio of particle size to interparticle spacing) are 
negligible, nor so small that the flow is choked. This will ensure that appreciable hydrodynamic flow takes place 
between the spheres, giving rise to the strain-dependent mobilities || that are used in (^). If the system parameters 
are as in [ p7| , the Reynolds number will be negligible, as required by our neglect of inertia, and the Peclet number 
large. Note that our model equations ([|) were formulated to describe the nature of distortions about a single 
crystalline domain. In particular, the instability towards clumping takes place only on large enough length scales. In 
a polycrystalline sample, if the size of the crystallites is too small, terms from the elastic energy in (||) could dominate 
instead. In addition, it is important that the sedimentation be steady, a requirement best met by working in the 
fluidised-bed geometry in which the particles constituting the crystal are on average at rest in the laboratory frame 
of reference, and the fluid flows vertically upwards past them. We would recommend starting with the suspension in 
the fully sedimented state, and then switching on the upward flow. Observations in |^8| suggest that strongly charge- 
stabilised crystalline suspensions appear stable whereas suspensions in a fluid state display the Crowley instability in 
a visible manner. We suspect that the instability is present even in the crystalline suspensions, but is masked either 
by finite crystallite size or by the logarithmically slow coarsening of domains. We predict that careful measurements 
of the time-evolution of the static structure factor, using particle-imaging or ultrasmall-angle scattering techniques, 
should reveal a weak large-scale modulation of the particle concentration, with characteristic wavevector normal to 
the sedimentation direction and decreasing logarithmically in time. 
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